Electron-electron interactions and charging effects in graphene quantum dots 
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We analyze charging effects in graphene quantum dots. Using a simple model, we show that, when 
the Fermi level is far from the neutrality point, charging effects lead to a shift in the electrostatic 
potential and the dot shows standard Coulomb blockade features. Near the neutrality point, surface 
states are partially occupied and the Coulomb interaction leads to a strongly correlated ground 
state which can be approximated by either a Wigner crystal or a Laughlin like wave function. The 
existence of strong correlations modify the transport properties which show non equilibrium effects, 
similar to those predicted for tunneling into other strongly correlated systems. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73.43.Lp 
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I. INTRODUCTION 

Graphene has attracted a great deal of attention, be- 
cause of its novel fundamental properties and its poten- 
tial applicationai. The interest on graphene devices has 
motivated recent research on the transport properties of 
small devices^i^i^i^. Features such as charging effects and 
quantum confinement are of crucial importance for their 
understanding^. The confinement of electrons, and the 
observation of Coulomb blockade effects has already been 
demonstrated experimentallj*i'2.. Note that the confine- 
ment of electrons in graphene is not trivial, due to the 
Klein's paradoo^, which makes potential barriers trans- 
parent for normally incident quasi-particles. Electrons 
in graphene can be confined, however, by exploiting the 
angular dependence of scattering at a barrier—. 

For graphene layers, electron-electron interaction is 
usually neglected, including works on localizationii even 
though disorder enhances the effect of interactionji^ The 
reasoning for this is to assume a "normal" ground-state 
at zero doping - characterized by a semi-metal. Be- 
cause the kinetic and interaction energy equally scale 
with the carrier density, the interaction does not be- 
come important at finite doping, either. It is thus well 
agreed on that at finite doping electron-electron inter- 
action can be treated within the random-phase approx- 
imation (RPA).— Nevertheless, at the Dirac point 
RPA seems to fail leading to a novel plasmon mode in 
graphene;i^ Also in a quantum dot, we find that electron- 
electron interaction has to be treated differently for dop- 
ing regimes close to and away from the Dirac point. 

The main part of this work is the prediction and 
characterization of strongly correlated few-electron states 
in graphene quantum dots. Similar studies have 
been performed previously for semiconducting quantum 
dotsi^iiiii^i^. In order to obtain strongly correlated 
ground states the Coulomb interaction has to dominate 
over the other energy scale, namely the shell structure of 
the single particle spectrum determined by the confine- 
ment. This is typically achieved by either using strong 
magnetic fieldsiiii^, so that the single particle levels form 



highly degenerate Landau levels, or using rather weak 
confinement."'^^ Interestingly, strongly correlated states 
naturally arise already in small graphene quantum dots 
even without magnetic field. The reason for that is the 
appearance of a highly degenerate zero energy band of 
surface states, which is strongly affected by Coulomb in- 
teraction. Close to half filling these states are occupied 
by few electrons which are strongly correlated and can 
be approximated by a Laughlin like wave function or al- 
ternatively by a quasi one-dimensional Wigner crystal. 

The paper is organized as follows. In section II, we 
present a simple model which allows us to describe qual- 
itatively the charging of a graphene dot. In section III, 
we show that the charging properties of the graphene dot 
are in agreement with Coulomb blockade theory when the 
Fermi energy is far from the neutrality point. Thereafter, 
we show in section IV that close to the neutrality point 
charging effects are strongly modified by the presence of 
midgap states, associated to the edges^^. We show that 
electrons occupying these midgap states form a strongly 
correlated state which is characterized in detail. In sec- 
tion V, we then discuss implications for transport prop- 
erties. We close with conclusions and outlook. 



II. THE MODEL 

The linearized tight-binding Hamiltonian for a 
graphene sheet with circular symmetry is given by 



VF 





'isdr 



isdr 




(1) 



where s = ± determines the valley. We assume that 
the dot is ballistic, i.e., with no internal disorder. The 
general solutions with energy efe = ztvpA: are of the type 



(2) 



with Jm{x) denoting the m-th Bessel function. The dot 
has a circular shape with radius R. The circular symme- 
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try of the dot allows us to classify the solutions according 
to their angular momenta. 

In order to analyze the possible role of surface states, 
we assume that the boundary conditions at the edges 
are those appropriate for a zig-zag graphene edge ending 
always on the same lattice site^S 



. 



(3) 



The boundary condition is not experimentally realizable 
for a circular dot, though it enables a detailed analysis of 
the interplay between Coulomb interaction and surface 
states. For the chosen boundary condition the wave vec- 
tor is quantized by fc = z„m/_R, where Znm denotes the 
n-th root of the m-th Bessel function, Jmiznm) = 0. In 
addition to the finite energy states given in Eq. ^ the 
boundary condition allows for surface states which can 
be written as 



*f(r,0) 







^^2(m + l) ' 



(4) 



with m > to guarantee normalizability. Note, that for 
the surface states the angular momentum is given by sto, 
see Eq. Q and that these functions have an analytical 
dependence on either z = x + iy oi z — x — iy. Discrete 
lattice effects impose a maximal (absolute) value on the 
angular momentum of order rrimax ~ R/a, where R is 
the radius of the disk and a is a length comparable to 
the lattice spacing. 

Charging effects arise from electron-electron interac- 
tion which is generally described by 



Ha = 



1 



n<n' 



(5) 



The total Hamiltonian is given by the sum of Eqs. ([T]) and 
(O. We note that both parts scale as 1/R. Furthermore 
in graphene (e^/47reoe)//iw_F ~ 1, so that both single par- 
ticle and interaction energies can be expressed in units of 
{e^ / iireoeR) — hvp/R, as will be done throughout this 
paper. 

Charging effects are mostly determined by the over- 
all geometry of the dot, so that the lack of disorder in 
the model described here does not change qualitatively 
the main features of Coulomb blockade. Graphene dots 
have, most likely, rough edges. Hence, the possible sur- 
face states are confined to certain regions of the edges. 
The model overestimates the number of surface states of 
a given dot. On the other hand, wave functions localized 
in the angular coordinate 9 can be built from the wave 
functions in Eq.® or Eq.Q. A dot where the edge has 
a region of size I of the zig-zag type has states localized 
at the edge with an angular width A6' ^ l/R. These 
states will be approximately described by superpositions 
of states with angular momenta m < R/l. Hence, when 
R/l » 1 and l/a ^ 1, these states, which will change 
over distances larger than the lattice spacing, will be well 
described by superpositions of the states derived from our 
continuum model. 



III. CHARGING EFFECTS AWAY FROM THE 
DIRAC ENERGY 

We analyze the effects induced by increasing the num- 
ber of electrons in the dot using the Hartree approxi- 
mation. The self-consistent Hartree potential describes, 
within a mean field approximation, the screening of 
charges within the dot. We assume that a half filled dot 
is neutral, as the ionic charge compensates the electronic 
charge in the filled valence band. Away from half filling, 
the dot is charged. Then, an electrostatic potential is in- 
duced in its interior, and there is an inhomogeneous dis- 
tribution of charge^i. We describe charged dots by fixing 
the chemical potential, and obtaining a self-consistent so- 
lution where all electronic states with lower energies are 
filled. The Hartree approximation should give a reason- 
able description when Coulomb blockade effects can be 
described as a rigid shift of the electrostatic potential 
within the dot- -'-?. 

The Hartree potential needs to be calculated self con- 
sistently which must be done numerically, despite of the 
simplicity of the model. The Dirac equation for each an- 
gular momentum channel is discretized, and an effective 
tight binding model is defined for each channel. Details 
are given in Appendix [X] 

The conservation of the angular momentum allows for 
the possibility of solving dots with a large number of 
electrons. Typical results for dots charged with elec- 
trons or holes away from the Dirac energy are shown 
in Fig. [TJ The calculation has been done in a discrete 
lattice with N — 100 sites (see the Appendix). The 
Hartree potential changes little within the dot and, to 
a first approximation, the deviation from neutrality of 
the dot can be approximated by a rigid shift of the elec- 
trostatic potcntiaPV 



IV. CHARGING EFFECTS NEAR THE 
NEUTRALITY POINT 

The Hartree calculations mentioned above fail to give a 
self consistent solution when the surface band is partially 
occupied and a more advanced treatment of the interac- 
tion has to be applied. This also implies that deviations 
from conventional Coulomb blockade can be expected in 
this regime. 

Instead of treating the interactions within a mean field 
approach we therefore employ the method of configura- 
tion interaction to fully take into account all correlations 
within the truncated Hilbert space of surface states. The 
truncation of the Hilbert space can be justified by the 
energy gap to extended states of finite energy which in 
our model is given by lAhvp I R- In principle, the effect 
of the extended states can be added to the following anal- 
ysis as a perturbation, but we do not expect qualitative 
changes of our main conclusions. 

In the following we deal, therefore, with a few-electron 
problem and consider the eigenspectrum of N interact- 
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FIG. 1: (Color online). Electronic structure of a quantum 
dot in the Hartree approximation. All energies are in units of 
hvpN/R with A'' = 100 and the position in units of r/R. 
Top part: Electron energies as function of the angular mo- 
mentum, for one valley in the Brillouin Zone, and different 
values of the chemical potential. Left: ep = 0.5. Center left: 
ep = 0.25. Center right: ep = -0.25. Right: ep = -0.5. The 
total number of states per valley is 12200. The number of oc- 
cupied states per valley is 6155 (left), 6122 (center left), 6066 
(center right), and 6038, right. Bottom part: Hartree poten- 
tial (top), and charge density (bottom) for the four values of 
the chemical potential considered in the top part: Black (solid 
line), ep = 0.5. Red (dashed line), ep = 0.25. Green (dash- 
dotted line), ep = —0.25. Blue (dotted line), ep = —0.5. 



ing electrons occupying surface states. The interaction 
is described in Eq. Since the screening of electron- 
electron interaction is known to be poor close to half 
filling, it seems sensible to consider a long-ranged inter- 
action rather than a (point-like) Hubbard interaction. In 
addition to the particle number, the few-electron wave 
function can be characterized by the valley-polarization 
Iz — (in absence of inter- valley scattering) , the 

total angular momentum M — J2n ^nmn as well as by 
the quantum number S'^,Sz for the total spin. In the 
following we limit ourselves to valley and spin polarized 
solutions. While spin polarized electrons cannot interact 



with each other via a point-like Hubbard interaction (due 
to Pauli principle), the long ranged Coulomb interaction 
will give rise to highly correlated spin-polarized states as 
shown in the rest of the paper. 

Electron-electron interaction tries to maximize the dis- 
tance between the electrons, which leads to a correlated 
ground state. The few-electrons ground state for Af ^ oo 
is given by a classical Wigner crystal where the N elec- 
trons are localized at r = i? and 6'„ = 27rn/iV, thus min- 
imizing the Coulomb energy, see insets in Fig. [51 It is 
important to note that due to the localization, the trun- 
cation of the Hilbert-space to include only surface states 
is still (in fact, even better) justified in the presence of 
electron-electron interactions. 

Surface states are characterized by only populating one 
of the two sub-lattices and thus avoiding the kinetic en- 
ergy due to nearest-neighbor hopping t. However, next- 
nearest neighbor hopping t' « t/10 « 0.3eV connects 
sites within the same sub-lattice so that surface states 
gain some finite kinetic energy and the zero energy band 
becomes dispersive. This kinetic term delocalizes the 
wave-function of the surface states and leads to a stable 
few-electron ground-state with finite angular momentum 
Mq. From ReffSi, the kinetic energy due to next-nearest 
neighbor hopping reads t'a^p^, with a the lattice spacing 
and p the momentum operator. As shown in appendix [B] 
the Hamiltonian for next nearest neighbor hopping Hkin 
can be written to lowest order perturbation in t' as 



hvp 3a 
R lOR 



^TO(m -I- l)cj„c„ 



This kinetic term competes with the Coulomb interac- 
tion, since it reduces the Coulomb correlations of the 
ground state. This competition is also visible in the de- 
pendence of the few-electron energy on the total angular 
momentum as shown in Fig. [2l In the absence of next- 
nearest neighbor hopping (solid line) the energy decreases 
with increasing angular momentum (except for oscilla- 
tions discussed below), since states of higher angular 
momentum have lower Coulomb energy. However, when 
next-nearest neighbor hopping is included then the occu- 
pation of states with large angular momentum is hindered 
and the energy as function of angular momentum shows 
a well defined minimum. We note that the ratio between 
kinetic energy and Coulomb energy increases with de- 
creasing dot size (the numerical calculations are done for 
R = 22nm). Consequently for smaller dots the angular 
momentum of the ground-state decreases. In the studied 
subspace of valley and spin polarized electrons, the mini- 
mal angular momentum is given by Mmin = N(N — l)/2. 



A. Trial functions for the correlated ground-state 

The lack of well converged Hartree solutions, which 
are given by Slater determinants, imply that the wave- 
function which describes the surface states in the pres- 
ence of charging effects is strongly correlated. We have 
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chosen two ansatze which are compared to the numeri- 
cahy exact solution. 



1. Laughlin wave function 

The appearance of a partially filled degenerate energy 
band separated by an energy gap from the rest of the 
spectrum strongly resembles Fractional Quantum Hall 
physics. However, now the zero energy band is caused by 
the boundary condition and the gap is due to the confine- 
ment rather than due to high magnetic fields. Not only 
the band structure is similar in both systems, but also 
the form of the one-particle states is similar. Both the 
surface states as well as the orbitals of the lowest Landau 
level (in symmetric gauge) depend on z™. 

This analogy can be used to propose a trial wave func- 
tion much like Laughlin's original wave function for the 
ground-state in the Fractional Quantum Hall regimo^S 



(6) 



with p odd to ensure antisymmetry and C a normalization 
constant. These wave functions have a well defined total 
angular momentum, M = pN{N — l)/2. For p — 1, this 
is the minimal possible angular momentum of A^-fuUy 
polarized electrons occupying surface states and the trial 
wave function (which in this case is given by a single 
Slater determinant) is the exact eigenstate. With in- 
creasing value oi p — 1, 3, 5, • • • the correlations increase 
and the wave function is given by an increasing number 
of superposed Slater determinants, much like Laughlin's 
original wave function for the Fractional Quantum Hall 
stated. The Laughlin like wavefunction in Eq. ^ is 
a parameter free trial wavefunction that conserves the 
present symmetries (i.e. total angular momentum) and 
that can be uniquely expressed in the subspace of sur- 
face states. Furthermore the factors {zi — Zjf create 
extended holes around each electron, which minimizes 
the Coulomb energy and explains the good agreement 
between the trial wave function and the numerically cal- 
culated ground state. We note however, that we use the 
similarity between the studied system and the Fractional 
Quantum Hall effect only to get a trial function for the 
ground state, while we do not analyze the similarities 
between both systems in the excitation spectrum. 



trial function as 



mi, m2, ■ ■ ■ , mjY 



(7) 



where the operator c]„ creates a state with momentum 
Tfin- Note that due to the constraint imposed by the 
antisymmetry requirement the wave function can only 
be defined for total angular momenta of the form M = 
N{N — l)/2 + jN where j is a positive integer. While 
the phase factor guarantees for the angular correlations 
in the wave function, we use the factor (rn„ -t- 1)™ to 
optimize radial correlations (note that the normalization 
constant of a surface state is proportional to ^/rr^ + l). 
For a strictly one dimensional system the usual definition 
of a quasi-classical Wigner crystal implies w = 0. In the 
numerical calculations, we choose w such that the wave 
function optimizes the ground-state energy [w ~ 2 for 
low M and t'). 



B. Correlations 

For AI —> oo and without next-nearest neighbor hop- 
ping, the quantum mechanical configuration approaches 
the classical one, which minimizes the Coulomb energy by 
pinning the electrons at r = i? and = 2T:n/N. How- 
ever, due to the rotational symmetry of our problem, the 
ground-state is a superposition of all orientations such 
that there is a constant density distribution around the 
circumference. To characterize a Wigner crystal or more 
generally a density correlated system, one thus has to 
look at the density-density correlation function 



Cf^(ro,r) 



{N,M;0\- 



,6{r- 



N{N- 1) 



-\N,M;0) 
(8) 



where \N,M;Q) is the ground state of the N particle 
system to fixed angular momentum M and i,j g 1, .., N. 



Effect of Disorder 



2. Wigner crystal 

An obvious alternative to the FQHE like wave func- 
tion described above is that of a Wigner crystal. The 
surface states are maximal at the border of the dot and 
the system resembles a one-dimensional system. In order 
to minimize the Coulomb energy, it is therefore favorable 
to superpose the wave functions in such a way that elec- 
trons are maximally separated in angle. We write such a 



We can use the same truncated basis to study disorder 
due to the roughness of the edges. Due to the flat disper- 
sion in the absence of disorder, single particle states tend 
to localize near imperfections of the edge, however one 
can show that the degeneracy of the zero energy states 
is only reduced by the number of impurities, which can 
be assumed to be much smaller than the number of sur- 
face states. The correlated state found in the presence 
of interactions will also be pinned by disorder, leading to 
glassy featuresS^. 
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FIG. 2: (color online). Ground-state energy (in units of 
for A*' = 2, 3, 4, 5 electrons occupying surface states 



as function of the total angular momentum M with (red solid 
line) and without (blue dotted line) next-nearest neighbor 
hopping t' (assuming R — 22nm). For t' = and M — > cxa 
the energy approaches the classical energy of point charges 
on a disc. The classical configuration is shown in the inset 
and the classical energy is indicated by the constant dashed 
line. 




FIG. 3: (color online). Ground-state energy for A'^ = 3 elec- 
trons with (full) and without (dotted) next-nearest neighbor 
hopping t' in comparison with the one obtained from the trail 
functions. The right hand side shows energy differences. En- 
ergies are in units of -r^ — 



D. Numerical results 

Fig. [2] shows the energy of the lowest lying spin- and 
valley-polarized eigenstate of each total angular momen- 
tum M for TV = 2,3,4,5 electrons occupying surface 
states. The energies are obtained by numerically di- 
agonalizing the few-particle Hamiltonian in this sub- 
space. The dotted line in Fig. [2] shows the results if 
next nearest neighbor interaction is neglected, so that 
the total Hamiltonian consists of the Coulomb interac- 
tion, only. We note two main features in that case. First, 
the energy oscillates as function of the angular momen- 
tum with local minima at M = Mmin{N) + jN, where 
Mmin{N) = N{N — l)/2 denotes the minimal angular 



momentum of N- spin and valley-polarized electrons and 
J is a positive integer. Only at these angular momenta 
the angular correlations between the electrons can be 
fully developed. This can be seen in the correlation func- 
tions discussed below and in the fact that only for these 
distinct angular momenta a Wigner trial function can be 
constructed. The second feature visible in Fig. [2] is that 
the energy generally decreases with increasing angular 
momentum for t' = Q and it finally reaches the classical 
limit corresponding to N point charges on the dot. The 
classical configurations are shown in the insets, and their 
energies are indicated by the constant dashed lines. 

For finite t' (see solid line in Fig. [2]), the Hamiltonian 
is supplemented by a kinetic term given in Eq. ([6]) that 
competes with the Coulomb interaction. Since the kinetic 
energy of surface states increases quadratically with their 
angular momentum the cost in kinetic energy exceeds for 
large total angular momenta the gain in Coulomb energy 
connected with an increase in angular momentum. Thus 
the A^-electron system now has a ground-state with well 
defined angular momentum Mq. The ratio between the 
kinetic term and the Coulomb energy grows for decreas- 
ing dot sizes, which also leads to a decrease in Mq. 

In Fig. Owe compare the numerically obtained energies 
for = 3 electrons with that of the two trial functions 
described above. The data for the Laughlin-like wave 
function (defined in Eq. ^) is indicated by squares while 
the data for a Wigner-crystal-like wave function (defined 
in Eq. ([7])) is labeled by filled circles. First, we note 
that the energies of both trial wave functions differ by 
less than 1% from the numerical data. As noted above 
the Wigner-crystal-like wave function can be constructed 
for each angular momentum where the few-electron en- 
ergy shows a minimum. In contrast a Laughlin-like wave 
function only exists for each TV — 1-th minimum. It is 
interesting to note that the Laughlin-like wave function 
becomes better for finite t' than for t' = 0. 

In Fig. [3] we optimized the free parameter w in the 
Wigner wave function for each M separately, which leads 
to this extremely good agreement with the exact data 
for both zero and finite t' . We note, however, that the 
optimal value was w k, 2 for all M in the case of = 0, 
while we strongly increased w with increasing M for finite 
t'. 

Fig. [3] shows the density plot of the exact, sym- 
metrized density-density correlation function C|^(r) — 
E'iLiC^jiR,i2n/N;r) for N = 2,3,4,5 electrons and 
for AI = N + Minin. The A^-fold symmetry, which is 
typical of a 1-d Wigner crystal is clearly seen. We note 
that also the trail wave function show these correlations, 
which explains the good agreement of its energies with 
the exact one. 

In Fig. [5l the angular correlations along the perimeter 
of the dot is shown for = 3. An electron is fixed 
at 6* = and r — R and the probability of finding 
another electron at a given angle is plotted. The left 
hand side of Fig. [5] illustrates that correlations are maxi- 
mally developed at the distinguished angular momentum 
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FIG. 4: (color online) Density plot of the exact, symmetrized 
density-density correlation function C'£(r) for A'^ — 2,3,4,5 
particles and total angular momentum M = N + Mmin- 



N = 3, f=0 N = 3.M=48 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

e ; In e ; In 



FIG. 5: (color online) Angular correlations C^f{R,0; R,9) for 
A'' = 3 along the perimeter of the dot for various total angu- 
lar momentum M (left hand side) and various next-nearest 
neighbor hopping t' (right hand side). 



M ~ jN+Mmin (here j — 3), while for other angular mo- 
menta the correlations are washed out. On the right hand 
side of Fig.[5l we see that the density-density correlations 
are more pronounced for higher angular momentum (here 
j = 15) while the kinetic energy t' reduces these correla- 
tions, which again is a manifestation of the competition 
between Coulomb interaction and next-nearest neighbor 
hopping. 

V. TRANSPORT PROPERTIES 

The addition of one electron to the dot, in the regime 
where the surface states are partially occupied, not only 
charges the dot and shifts the electrostatic potential, but 
changes the correlated wave function as well. Hence, one 
expects a correction to the local density of states in the 
dot, which is energy dependent, in a similar way to An- 
derson's orthogonality catastrophe^!, or the singularity 
in the X-ray core level photoemissio n^^i^^ . Such Fermi 
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FIG. 6: (color online). The (unnormalized) spectral function 
An{uj) for = 3 (left hand side) and N = 4 (right hand 
side) with respect to the ground-state energies E^' and -Bq'^^, 
respectively. The full, dashed, and dotted lines correspond to 
different maximal angular momentum. 



edge singularities have also been discussed in relation to 
transport in quantum dots and nanotubes^SiSiiHi^^ ^ 

The correlated state which describes the surface states 
of the graphene quantum dot resembles a one dimen- 
sional system localized along the surface. In this respect, 
the tunneling into this state can also be analyzed within 
the related framework of tunneling into correlated one 
dimensional metals"^"*. In this case, and in those describe 
before, one expects that the tunneling density of states 
of the dot will be described by a power law. We have 
computed numerically the spectral function 

"Imax 2 

^|(iV - l,Afo;0| ^ c„|iV,M;n)| 

A/,n m—0 

X S{e!!^'' - eI^-''''- - Lo) (9) 

where Cm annihilates a particle with angular momentum 
TO. Next nearest neighbor hopping t' causes a finite total 
angular momentum Mq of the N— 1-electron ground state 
and due to momentum conservation the angular momen- 
tum of the N-electron state is given by M = Mq -I- to. 
We restrict the m-summation by an upper angular mo- 
mentum. 

Results for the spectral function are shown in Fig. [5] 
which are characterized by a sharp peak, reminiscent to 
the delta peak of the non-interacting system. Due to 
the electron-electron interaction, this peak is smeared out 
and decays as a power law decay, in qualitative agreement 
with the arguments mentioned above. There is a clear 
convergence for low energies as function of the maximal 
angular momentum. 

VI. SUMMARY AND OUTLOOK 

We have presented a simple model of a graphene quan- 
tum dot, suitable for the analysis of interaction effects. 
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We show that Coulomb blockade effects are similar to 
those in other systems when the chemical potential is far 
from the neutrality point. 

The Dirac equation which describes the electronic 
states of graphene allows for the existence of midgap 
states, near defects or surfaces. The presence of these 
states changes qualitatively the properties of the dot in 
the Coulomb blockade regime. As the kinetic energy of 
these states is nearly zero, the resulting wave function is 
mostly determined by the interaction, and deviates sig- 
nificantly from a single Slater determinant. In order to 
describe correlations beyond mean field we employed the 
method of configuration interaction within the subspace 
of surface states. Since it is known that screening is weak 
in the described case close to half filling, we considered 
the electrons to interact with each other via the long- 
ranged Coulomb interaction in contrast to a point-like 
Hubbard interaction studied for example in Rcf.'^-. 

Making use of the simple analytical form of the surface 
states, we have identified two possible correlated wave 
functions which are in good agreement with few-particle 
exact calculations: a wave function similar to that pro- 
posed by Laughlin for the Fractional Quantum Hall ef- 
fect, and another describing a Wigner crystal. These 
results indicate the existence of strong correlations, al- 
though they do not allow us to analyze the existence of 
an incompressible electron liquid in the thermodynamic 
limit. We note that the correlations present in the spin 
polarized states studied here, arise only for long-ranged 
interactions, while they are absent if an effective point- 
like Hubbard interaction is considered. We expect the 
few particle states to be pinned by disorder at the edges 
(note, however, that their extension can be comparable 
to the dot size, so that the pinning will not be large). 

The transport properties in the regime where the 
midgap states are partially occupied will deviate from 
that observed in other quantum dots. The strongly corre- 



lated nature of the wave function implies that non shake 
up effects will suppress the tunneling density of states. 

An interesting extension of this work concerns the 
valley- and spin-degree of freedom which will be ad- 
dressed in a subsequent publication. 
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APPENDIX A: DISCRETIZATION OF THE 
DIRAC EQUATION 

The Dirac equation for angular momentum / can be 
written as two coupled one dimensional differential equa- 
tions: 

V{r)'il;A{r) + vp (^idr ± ipsir) = e^^Air) 

VF (^idr T ipAir) + V{r)'ipBir) = etpsir) 

(Al) 

where the two signs correspond to the two Dirac points. 
We now analyze a given Dirac equation. Extension to 
the other valley is straightforward. Equation (jAl|) can 
be written as: 



V{r) li^Air) + i^sir)] + vp [ idr + [V'a(»') + ipBir)] - ivp J" ^ [tpAir) - i^sir)] = e [ipA{r) + V's(^)] 

Ir J Ir 

V{r) {^A{r) - ^B{r)\ - «f ( ^dr + ) [VM(r) - ^^(r-)] + «x'f^^ \^A(r) + ^B{r)\ = e {^A(r) - ^B{r)\ (A2) 

' 2r / 2r 



I 

We define and we obtain: 



Vii(r) 



i>A{r) - i>B{r) 



(A3) 



+ \ ~ 

V(r)tpi{r) +ivFdrtpi{r) - ivF— ip2{r) = etpi{r) 

2r 

21 + 1 ~ 

V(r)tp2{r) - ivFdrtp2{r) + ivp—z ipiir) = etp2{r) 

Ir 



(A4) 
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l+(21+lV(4n) 



l-(21+l)/(4n) 



a|i+2 b|,+2 



and: 



2tt 



dO- 



(A7) 



N 



In terms of the original Dirac equation, the energies are 
expressed in units of hv^ / R and the parameter wq is given 
by Wo = {e^ l'^T:eot)/hvF « 1. 



FIG. 7: (color online) Top: Sketch of the discretization of the 
Dirac equation in radial coordinates used in the text. Bot- 
tom: doubled chain used in the calculations, in order to avoid 
spurious effects at n = 1. 



A set of discrete equations which, taking the contin- 
uum limit, lead to Eq. (IA4p is: 



21 + 1 

An 
21 + 1 

in 



1 



21 + 1 



An 



4n ' " 



K-i - 1 - 



+ Vna„ = ea: 



(A5) 



This set of equations is formally equivalent to a dimer- 
ized tight binding chain, as schematically shown in Fig.[71 
These chains admit zero energy estates localized at the 
ends, when the last hopping is smaller than the previous 
one. In order to avoid the formation of a spurious level at 
the center of the dot, n = 1, the chain is doubled, as also 
shown in Fig. [71 The Coulomb potential is discretized as: 



N 



m— 1 



(A6) 



APPENDIX B: KINETIC ENERGY DUE TO 
NEXT-NEAREST NEIGHBOR HOPPING 

Due to next-nearest neighbor hopping t' ~ O.lt, the 
initially flat band of surface states becomes dispersive. 
From Ref.— , the kinetic term due to t' is given by T = 
jt'a^p^ —jt'a^A. As for the Coulomb interaction we 
restrict our Hilbert space to the surface states V'm {r, 9) — 
"^f^+^r.e) defined in Eq. dH): 



(mlrln) = 



d rip:^Aipm ^ - J d rV?/',*„V'(/'m + 27r [■il^^rdripm]r=B^ 

In the second row we used partial integration leading 
to the boundary term (second term on right hand side). 
Including next nearest neighbor hopping this boundary 
term has to vanish, while the general form of the wave- 
function is assumed to change only close to the boundary. 
We thus only keep the first term which results in 



(mlTln) = S„ 



9aH' 



hvp 3a 



m(m + 1) 



^ A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

2 B. Huard, J. A. Sulpizio, N. Stander, K. Todd, B. Yang, 
and D. Goldhaber-Gordon, Phys. Rev. Lett. 98, 236803 
(2007). 

^ M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 

Rev. Lett. 98, 206805 (2007). 
* J. R. Williams, L. DiCarlo, and C. M. Marcus, Science 

(2007). 

^ J. B. Oostinga, H. B. Heersche, X. Liu, A. F. Morpurgo, 

and L. M. K. Vandersypen (2007), arXiv:0707.2487. 
® F. Munoz-Rojas, D. Jacob, J. Fernandez-Rossier, and J. J. 

Palacios, Phys. Rev. B 74, 195417 (2006). 

F. Sols, F. Guinea, and A. H. Castro Neto (2007), 

arXiv:0705.3805. 
® J. S. Bunch, Y. Yaish, M. Brink, K. Bolotin, and P. L. 

McEuen, Nano Lett. 5, 2887 (2005). 
^ M. L Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 

Physics 2, 620 (2006). 
° P. G. Silvestrov and K. B. Efetov, Phys. Rev. Lett. 98, 

016802 (2007). 



" L. Aleiner and K. B. Efetov, Phys. Rev. Lett. 97, 236801 

(2006) . 

T. Stauber, F. Guinea, and M. A. H. Vozmediano, Phys. 
Rev. B 71, 041406(R) (2005). 
" B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New J. 
Phys. 8, 318 (2006). 

E. H. Hwang and S. DasSarma, Phys. Rev. B 75, 205418 

(2007) . 

S. Gangadharaiah, A. M. Farid, and E. G. Mishchenko 
(2007), arXiv:0710.0622. 
^® A. Ghosal, A. D. Guclu, C. J. Umrigar, D. Ullmo, and 
H. U. Baranger, Phys. Rev. B 76, 085341 (2007). 
G. Yannouleas and U. Landman, Phys. Rev. B 66, 115315 
(2002). 

G. Yannouleas and U. Landman, Phys. Rev. B 70, 235319 
(2004). 

S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

Note, however, that electrostatic considerations in two 



9 



dimensions lead to a charge distribution which varies 
smoothly over the dot, but has a p(r) cx R — r di- 
vergence at the edges. 

D. V. Averin and K. K. Likharev, in Mesoscopic Phenom- 
ena in Solids, edited by B. L. Altshuler, P. A. Lee, and 
R. A. Webb (Elsevier, Amsterdam, 1991). 
H. Grabert and M. H. Devoret, eds.. Single Electron Tun- 
neling (Plenum, New York, 1992). 

N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 
Rev. B 73, 125411 (2006). 

R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

A. L. Efros and B. L Shklovskii, J. Phys. C: Sol. St. Phys. 

8, L49 (1975). 

P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 

P. Nozieres and C. T. de Dominicis, Phys. Rev. B 178, 



1097 (1969). 

G. D. Mahaii, Many Body Physics (Plenum Press, New 
York, 1993). 

M. Ueda and F. Guinea, Zeis, fiir Phys. 85, 413 (1991). 

E. Bascones, C. P. Herrero, F. Guinea, and G. Schon, Phys. 
Rev. B 61, 16778 (2000). 

D. A. Abanin and L. S. Levitov, Phys. Rev. Lett. 93, 

126802 (2004). 

F. Guinea, Phys. Rev. Lett. 94, 116804 (2005). 

C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992). 

J. Fernandez-Rossier and J. Palacios, Phys. Rev. Lett. 99, 
177204 (2007). 



